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The ground-state electron density of a one-dimensional spin-orbit coupled quantum dot with a 
Zeeman term and strong electron interaction is studied at the fractional helical liquid points. We 
show that at fractional filling factors v — (2n -|- 1)“^ (with n a non-negative integer) the density 
oscillates with Nq/iz peaks. For n > 1 a number of peaks larger than the number of electrons No 
suggests that a crystal of fractional quasi-particles with charge ve (with e the electron charge) 
occurs. The reported effect is amenable of verification via transport measurements in charged 
AFM-coupled dot. 
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I. INTRODUCTION 


Quasi-helical electrons occur in spin-orbit coupled 
quantum wires when a magnetic field perpendicular to 
the spin-orbit field is appliecPH^. The magnetic field 
opens a gap at the degeneracy point of the spectrum, 
when kp = fcsoi with kp the Fermi wavevector and fcso 
the wavevector associated with the spin-orbit interaction. 
For not too strong magnetic induction, the system is es¬ 
sentially helical with right (left) moving electrons having 
spin down (up). 

These systems, object of recent intense theoretical in¬ 
vestigations, possess several unusual properties such as 
the occurrence of Majorana bound states at their edges 
when proximized with a superconductoi!^ or peculiar 
spin oscillations in the presence of magnetic impuritie^. 
Most striking are the effects induced by repulsive elec¬ 
tron interactions: as an example, the gap induced by the 
magnetic field is strongly enha nced by interactions and 
anisotropic spin properties occuJiESIIil 
Strong interactions generate even more interesting effects 
in quasi-helical wires, when the Fermi level is tuned below 
the degeneracy point. A model originally introduced to 
study fra ctiona l gapped phases in arrays of parallel quan¬ 
tum wireJi^IISl treated as Luttinger liquidJlsHIIl ]jag been 
recently applied to the four-channels Luttinger liquid de¬ 
scribing spin-orbit coupled quantum wirealSl. The com¬ 
bined effects of the magnetic field and electron-electron 
scattering are captured by a non-linear Zeeman coupling 
dressed by interactionJi^. Critical positions of the dot 
Fermi level are very relevant for fermions, characterized 
by a fractional filling factor with odd denominator 


kp 

kso 


1 

2n-\-l 


(n e N*). 


Here, interaction becomes resonant and for sufficiently 
strong interactio ns it is relevant in the sense of the renor¬ 
malization groujP^El! As a consequence, additional 


gaps open and fractional phases stabilize in the wire. 
They are accompanied by several signatures, such as 
the o ccurre nce of low-energy excitations with fractional 
charg^lSliSl which induce a non-integer quantization 
of the conductanc^isHini Also the shot noise becomes 
fractionalized^ and unusual properties of the chiral 
currents have been reported very recentlji^. 

What happens to ground-state properties such as the 
electron density of a quasi-helical quantum dot when a 
fractional phase occurs? This is a non-trivial question, 
since it can be expected that such a peculiar state may 
exhibit unusual fea tures, which wou ld comp ete with the 
conventional Friedel2212il and Wignei^^^lSMini oscillations. 
In any finite-size one-dimensional system, such as a 
quantum dot embedded into a quantum wire, electron 
density oscillations occur. A competition exists between 
oscillations with 2kp wavevector (Friedel oscillations), 
due to reflections at the dot edges, and oscillations 
with wavevector Akp (Wigner oscillations) due to strong 
electron repulsiorP^. The prevalence of either of the 
two is in general controlled by the ratio between the 
average electron repulsion and kinetic energy. In a dot 
with one charge and one spin degree of freedom, Friedel 
oscillations have No/2 peaks (with Nq the number of 
electrons, assumed even for simplicity) while Wigner 
oscillations display Nq peaks. In quasi-helical systems, 
Friedel oscillations in the main gap (an integer phase 
with ^ = 1) are predicted to display Nq peaks due to 
the locking between spin and momentum which halves 
the number of degrees of freedorrPH. 

In this paper, we address this point evaluating the elec¬ 
tron density of a quasi-helical quantum dot with open 
boundary conditions in a fractional phase u = (2n-|-l)~^. 
Along with the standard Friedel and Wigner terms and 
their higher harmonic^^, we consider contributions 
originating from the dressed Zeeman interactions. We 
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then evaluate the ground-state average electron density. 

Our main result is that in the fractional phase 
V = (2n -I- 1)“^ the most relevant oscillations of the 
electron density have wavevector 4(2n-|- l)fcF and exhibit 
(2n -I- l)iVo peaks. This is in striking contrast with 
the standard situation where for strong interactions A^o 
peaks occur, namely one for each electron. The presence 
of a larger number of peaks than the number of electrons 
suggests to interpret this result as the manifestation, 
in the ground state, of a crystal of fractional quasi¬ 
particles with each density peak corresponding to a 
fractional charge ve, with e the electron charge. Charge 
density oscillations can be experimentally detected, for 
instance in a linear transport experiment pe rformed 
while scanning the dot with a charged AFM tipPsEHM]^ 
We estimate the conditions to observe this effect in 
state-of-the-art sample^^lHlIl a measurement of such 
fractional density oscillations would then constitute a 
probe of the fractional states. 

The paper is organized as follows. In Sec.|TT]we introduce 
the Luttinger model of the dot with spin-orbit coupling 
and magnetic field. Focusing on odd-denominator reso¬ 
nances, we derive an effective Hamiltonian in the strong 
interactions regime and discuss the electron density op¬ 
erator. In Sec. |T^ our results are reported and discussed 
in details. Section Hvl contains our conclusions. 


II. THE MODEL 


fields ifnix), satisfy canonical commutation rules 

{4i^{x),dxi9^i{x')] = iTrSfj,^fj,rS{x — x') and can be repre¬ 
sented as 


,{x) = iy/% ^ 


r. — TTUaCt 


/2i 


nQ>0 




— 7rn„a 


1 X-V 


jll 


cos 


n„>0 


TTUqX 


ft,rig ' j I 


where a is a short-length cutoff and are canonical 
bosonic operators. The fermionic operator is expressed 
in terms of left and right components 

4',(x) = e-^^^^L,{x)] , (3) 


with s = ± the z direction of the electron spin. Here, 
fcp = 'xNq/{2£) the Fermi wavevector, with Nq the refer¬ 
ence number of electrons in the Fermi sea, assumed even 
for simplicity. The left- and right-moving fields Ls(x), 
Rs{x) satisfy Ls{x) = —Rs{—x) with 

Rs{x) = -, (4) 

V27ra 

Fs a Klein factor, the number of excess electrons in 
the s = ± spin sector with respect to the Fermi sea and 


<i)s(x) = 6p{x) F s9„{x) - 4>p{x) - S(j)cr{x) ■ 
Note that [$s(x), $s/(—x)] = 8i5s,s'f{x) with 


fix) 


1 

- arctan 
2 


sin 




g7ra/f _ gQg 
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27rx \ 
£ } 


(5) 

( 6 ) 


A. Zero magnetic field 

We consider a one-dimensional quantum dot of length 
i with open boundaries at x = 0 and x = 1. The dot 
is subject to a spin-orbit field acting along the negative 
d irec tion of the z axis. The single-particle Hamiltonian 
ijm 


due to the finite size of the dofP^. 

When no excess charges are considered, Ns = 0, the 
Fermi energy crosses the energy spectrurrP at the four 
wavevectors ±(fcF ± kso)- 

B. Magnetic field 


Hsp = + igsodx(T^ , (1) 

where m* is the effective electron mass, gso the strength 
of the spin-orbit coupling, with associated wavevector 
kso = w*5so, and 7, cr^ the identity and z Pauli matrix 
respectively (here and henceforth h = 1). 

Including forward electron-electron interactions, the 
Hamiltonian around the Fermi energy can be described 
as a two-channels Luttinger liquid 


A magnetic field along the x axis induces a Zeeman 
coupling V = ^p,Bg*B cTx = with B the field inten¬ 
sity, fj,B the Bohr magneton and (Tx the x Pauli matrix. 
In terms of the fermionic operators 



'I'+(x)4>_(x) -I- 4'(l(x)'I'+(x) 


(7) 


When 'I's(x) is expressed via Ls{x) and i?s(x) eight terms 
arise 


— idx(l>pf + 9ii [dxdp)^ 

9n 

( 2 ) 

with p, = p (a) the charge (spin) modes and Vp their 
group velocity. The parameter g^ controls interactions 
and Vp = vp/gp in the p sector, with vp the Fermi 
velocity. For repulsive interactions one has 0 < Pp = 
g < 1, while g^ = ^ when SU(2) invariance holds. The 




( 8 ) 


In the presence of electron interactions, each term Vjjx) 
is also dressed by backscattering interaction vert ice 

, -in __ 

Unix)= e2*'=''“L^±(x)i?±(x) aPl 


Vjix) ^ h/”)(x) = 17„(x)y,(x)C/„(x) , (9) 
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with corresponding Hamiltonian 


condition - see also Sec. mil 


„>0 j = l 

In bosonic form, 

cos 2g|”^a; — V20l^\x) 


( 10 ) 


( 11 ) 


where are interaction amplitudes cx B. The 


3 

o, (») 

(^Ti. ; 

Q 

^.7 

1 

( 7 „ - l)0p(x) -b 9a(x) - 4>a{x) 

(7n — Ijfep + kso 

2 

( 7 „ - l)<t>plx) -b 9a{x) -b falx) 

(jn — l)kF + kso 

3 

( 7 „ - 2)(l)p{x) -b 9„{x) 

(jn — 2)kF + kso 

4 

Xnfpix) -b 9cr(x) 

XnkF + kso 

5 

( 7 „ - l)4>p{x) - 9a{x) -b <t>a{x) - F^^Jx) 

(jn — Ijfcp — kso 

6 

( 7 „ - l)<j)p{x) - 9a(x) - (f)a{x) - FI^"\x) 

(jn — l)fcp — kso 

7 

Xnfpix) - 9a{x) - F^"\x) 

7nfcp — kso 

8 

(7n - 2)<?ip(a;) - 9„{x) - F^^^Jx) 

(jn — 2)kF — kso 


TABLE I. Bosonized expression of the eight Zeeman terms, 
dressed by electron-electron interactions and the correspond¬ 
ing wavevector. Here, y„ = 2n + 1 and Fq"'\x) = 2%/2(l — 

'yn)f{x). 


operators (x) and wavevectors are presented in 
Tab. 1^ where 7 „ = 2n -I- 1. We have omitted zero modes 
describing excess electrons beyond the Fermi sea, which 
are irrelevant for the forthcoming discussions focusing 
on the behaviour at the Fermi surface with Nq electrons. 


The first four terms 1 < j < 4 oscillate with a wavelength 
shorter than i and therefore give negligible contributions 
to th e Ham iltonian in the sense of the renormalization 
groupIMISl. 

On the other hand by suitably tuning fcso and kp terms 
with 5 < j < 8 can become spatially non-oscillating 
and thus resonant. In particular, if kp = fcso/(2n) 
(n > 1), the terms V^"'\x) and Vq'^\x) become reso¬ 
nant, while for kp = fcso/(2n -1-1) (n > 0), the terms 
Vy (x) and (a;) resonate. Terms with j = 5,6 

(j = 7,8) are predi cted to give rise to fractional states 
in the dotP^ESElIill with filling factor v = 1/ (2n) (zz = 
l/(2n+l)). 

The above resonance conditions can be rewritten in terms 
of the average density no = No/£ as 


no 


2m*gso 


V 


( 12 ) 


where we have re-inserted h for clarity. For a resonance 
with a given v, the right-hand side of Eq. (121 is fixed 
by the material parameters. By tuning the average 
particle density via No and i the dot can be brought into 
resonance. Note that, due to the finite dot length, kp is 
quantized so that for a given n and number of electrons 
No only certain values of fcso satisfy the resonance 


In the following, we will only consider odd denominator 
resonances, i.e. those at kp = fcso/(2n-|- 1). Once res¬ 
onating, the terms V’j'^\x) and b ecome rele¬ 

vant in the spirit of the renormalization groupP^^ when 
g Qc with gc = 3/(2n -|- 1)^. Thus, for sufficiently 
strong interactions one can retain only the dominant 
terms of Eq. (fTol) 


1 


V K. — f dx 
27r Jo 


Vl"\x) + vt+^\x) 

A N 

= ^ y dx cos[4/(x)]cos[277+(x)] , 


where A^"^ = = A and 


(13) 


V+ix) = [infipix) - e„{x)\ -b 27 „/(x) . (14) 


The Hamiltonian H — Ho + V can be decoupled into a 
massless and a massive term, H pe ho -\- Hm, with 



— +gr, {d^X-f 

I9ri 

— {d^h+f +9v {dxX+f 
.9v 

cos [4/(x)] cos [2?7 _|_(x)] . 


(15) 


(16) 


With the massless field r]-{x) = [yn4>p{x) + 6fj{x)]/y/2 
and canonically conjugated X±{x) = 

[ln^(^p{x) T (t>cr{x)] /V2. Further, we have 


Vp = 


1 + 7^ / 1+7^ 

< 1 9ri — ln\N 22 

4 5r, V 1 + 7n5^ 


9- 


(17) 


Note that identifying H pe ho -\- hM neglects small 
corrections cx {8^^+) (dxV-) and oc {dxX+) i^xX-) 
which renormaliz^ Vjj and grj. The massless sector is 
a standard Luttinger liquid with renormalized param¬ 
eter while the massive term is a sine-Gordon model. 


The dynamics of the massless sector is standard. Con¬ 
cerning the massive sector with a sine-Gordon term, the 
equation of motion for ri+{x) is 

—dlg+ix) + Acos[4:f{x)]sm[2g+{x)] =0. (18) 

9v 

In the scaling limit 5 «C 5c which we consider, A di- 
verges under renormalization group transformational^*!^ 
and the semi-classical solutions for the equation of mo¬ 
tion ’’pin” the field r]+{x) to the minima of the cosine 
term, i.e. 


'n+{x) = ^ + kTr (keZ), (19) 
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with the conjugated field x+{^) & strongly fluctuating 
one: (x^(a;)) —)■ oo where (...) represents the quantum 
average over the groun d state. This is the infinite mass 
limit of the modePsHH. 


When the field rj+{x) is pinned as in Eq. (191 a gap in 
the spectrum opens. For n = 0 (^ = 1), Eq. (131 is 
relevant for g < 3 and corresponds to the opening of a 
gap near the degeneracy point of the spin-polarized wire 
subbands. This constitutes the main gap of the theory, 
which requires no repulsive interaction among electrons 
to be formed, where a quasi-helical electron liquid devel- 
opJIHll_ 

More intriguing are the properties of the fractional gaps 
occurring for n > 1 at fractional fillings 


1 

2n -I- 1 


These gaps have no counterpart in a noninteracting 
theory and will be the ones addressed in the rest of the 
paper. 


C. Electron density operator 


In order to study the properties of the ground state 
density, we start by observing th at the den sity operator 
p{x) is composed of several term^^nHSEUsslET] 

p{x) = pim{x) + XfPf{x) + XwPw{x) + Xzpz{x) , (20) 

with the weights Xi {i G {F, W, Z}) free parameters to 
be suitably determined. 

Here, 

Plw(x) = ^ + -d^(f>p(x) (21) 

ITT 

is the long-wave part of the density. 

The contribution of the Friedel oscillationJSMUI] jg 

OO 

Pf{x) = p^jr\x) (22) 

m—l 


with 


P^\x) = cos ^2m[kFX — f{x)] — V2m(j)p{x)'^ ■ 


cos 


V2m(j)^{x) 

the Friedel harmonics at wavevector 2mkF (m > 1). 
Interactions induce Wigner oscillation^^^ISlHII] 


(23) 


Pw(a;) = ^ Pw^(a:) (24) 

m—l 


with harmonics 

p^\x) = cos ^Am[kFX — f{x)] — 2'/2m(j)p{x)'^ (25) 

at wavevector AmkFX. 


Finally, a contribution induced by all terms in Eq. (10) 
which do not become resonant 


OO 8 


P'Z.ix) = X! X! “ {Sm,nSj,7 + <5m.n+l<Ji,8)] P^^j , 

m—l j — 1 


with 


(26) 


Pzj^(a;)=cos 2q^J^\ -'/20^^\x) . (27) 


is considered. 

These latter terms are motivated by observing that, to 
the lowest order, an interaction term of the form vj"^\x) 
induces a perturbation Sp^ j (x) to the long-wave terrrf^ 




Ee — Eq 


where \£) represents a many-body excitation of the un¬ 
perturbed dot with ground state |0) and E£ is the cor¬ 
responding unperturbed energy. Note that terms with 
j = 7, m = n and j = 8, m = n + 1 must not be 
considered here, since they are explicitly included in the 
Hamiltonian. It can be shown that 




COS 


(m) 

q) >X 


(28) 


namely vj^\x) produces oscillations with wavevector 

qj™). The inclusion of the terms p^^j{x) in the density 
operator captures the effects of residual interactio ns not 
included in the effective Hamiltonian of the systenP^E^. 
The weights {i G {F, W, Z\) in Eq. ( [^ are free pa¬ 
rameters to be determined imposing suitable constraints 
on /3(x), namely open boundary conditions and the con¬ 
servation of the number of electrons. 


III. RESULTS 

A. Density oscillations and fractional 
quasi-part ides 


As discussed in Appendix]^ in the infinite mass limit 
when r]+{x) is pinned - see Eq. ( |19[ ) - the only terms con¬ 
tributing to the average density^ (a;)) have wavevector 
I ■ (with I G Z*) and scale as [K{x)f with 


K{x) = 


sinh (^) 


v/sinh2(if)+sin2(^) 


(29) 
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see Eq. (A17). When g ^ 0 (and thus —>■ 0) the most 
relevant terms can be selected and the density, takes the 
final form in Eq. (A23) 


{ .gr? ^ N 

1 - ^ cos [A-f^kpx - F,(x)] I , 

(30) 

where Fi{x) = F 2 {x) = F 3 {x) = 0, F 4 {x) = 

4(7n - and F 5 (x) = 4(7„ + l)/(x). The 

terms E 4 ^ 5 (x) oc f{x) are a direct consequence of the 
dot finite siz^^ and stem from the terms p^^j{x) and 
(x) respectively. 


In the main gap n = 0 (with 7 „ = 1), the only non¬ 
vanishing contributions t o p( x) are the lowest Wigner 
harmonic (x) - see Eq. (25 1 - and four terms stemming 
from pz{x). They all have wavevector dfcp. The density 
displays Nq oscillations corresponding to the expected 
Wigner oscillations of Aq electrons in a bojP^. 

For n > 1, however, the density exhibits (2n -|- l)Ao 
peaks. This is due to the contribution of the relevant 
non-zero terms in pz{x) and of the (2n -|- l)-th Wigner 
harmonic. 

Figure and Fig. show (p(x)) for respectively Nq = & 
and Nq = 10 particles, in the case of the lowest 
resonances n = 0, n = 1, and n = 2. As can be 
seen, the density shows the expected number of peaks, 
Np = {2n -I- l)Ao. Superimposed to the peaks is a 
slowly-varying modulation induced by the terms F 4 5 (x) 
in Eq. (30 1 . The envelope due to [K{xyf''^‘^ is also 
present, but its modulation effect is negligible in the 
strong interaction regime presented here. 


For n > 1, a density displaying Np > Nq peaks is a 
fascinating situation where the many-body wavefunc- 
tion of the Nq electrons splits up. This splitting can 
be suggestively interpreted as the crystallization of 
(2n -I- l)iVo fractional quasi-particles, with fractional 
charge e/(2n -1-1), which tend to distribute evenly in 
space in order to minimize the strong repulsive interac¬ 
tions giving rise to Np distinct peaks in the density. 

It is very interesting to observe that, although several 
predictions have been already made concerning fractional 
excitations in strongly interacting spin- orbit coup led 
wires in the presence of a magnetic fielcP^E^IIIIlII^ jjj 
this work we predict that a bulk ground-state property - 
the density p(x) - exhibits signatures of these fractional 
quasi-particles. 


In this context, very recently a similar phenomenon 
has been reported in which the density of strongly 
interacting edge states of a two-dimensional topological 
insulator with two-particle backscattering exhibits 2Nq 
peaks for Nq electrons, and has been attributed to the 
formation of a ’’fractional Wigner crystal”!^. 

We also expect that signatures of fractional quasi- 



0 X i 



FIG. 1. Plot of the p{x) (units £~^) for a dot with No = 6 
electrons and g = 0.1 at the resonance (a) n = 0; (b) n = 1; 
(c) n = 2. In all panels, a = 


particles show up in higher-order functions, such as in 
the density-density correlators^ (not shown). 


We close this section observing that in deriving Eq. (301, 
all the non-zero dressed Zeeman terms p^^j (x) have been 
weighted with the same factor Xz- However, since all 
the relevant terms have the same wavevector, we expect 
that the number of oscillations of the density is robust, 
with only a possibly different modulation pattern of the 
peaks due to different weights of terms with Fi{x) and 
Fq{x). It may be also speculated that within the frac¬ 
tional phase the high-harmonic Wigner oscillations have 
a weaker amplitude that the term pz{x). Also in this 
case nonetheless we conclude, in analogy to what dis¬ 
cussed previously, that our results are robust barring a 
slightly different modulation pattern of the peaks. 
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FIG. 2. Same as in Fig. but for No = 10. All other param¬ 
eters as in Fig. 


B. Discussion 


We comment here on the observability and stability 
of the predicted effect. 

Assuming InAs as a host materiaP^^^, tuning the 
resonance at ly = 1/3 for a dot with Nq = 10 electrons 
one requires a lengttP^ £ « 7.5 /im, while £ « 4.5 ^m 
is required for Nq = 6 - see Eq. As put forward by 
several authors, density oscillations in a quantum dot 
can be probed by scanning a charged AFM tip along 
the wire while performing a linear tra nsport exper iment 
with lateral source and drain contact^^lMIllIllIlIl, xhe 
tip induces a chemical potential shift proportional to the 
dot density, and a modulation of the co nductance which 
is connected to the density oscillation^^UllIllIlIl^ Xhe 
lateral resolution can be approximately estimated in the 


range oP^ESI Ki 50 nm and the average wavelength of the 
density oscillations is « 250 nm. Thus, such an experi¬ 
ment may detect the Np peaks of the density oscillations. 


The terms V^j{x) and can become relevant 

as long as their wavelength A = 2-Kq^'^ is larger than the 
dot length £, or 


I . 27r 

|7„/cf - Ksol < ■ 


(31) 


Thus, the number of dot electrons and its length must 
be calibrated in order to satisfy the above condition. 
Note that Eq. (311 may allow for some flexibilty for the 


resonance with v = 1/3: supposing for instance that 
— kso = 0 for Nq dot electrons, one finds that 
the condition in Eq. (31) is still satisfied for Nq zL 1 


electrons in the dot. This fact does not hold however for 
resonances with n > 2, for which then the tuning of dot 
parameters is more critical. 

The arguments developed in this paper rely on the infi¬ 
nite mass approximation for the sine-Gordon problem, 
when the gap becomes the largest energy scale. Finite 
size effects limit the flow of the renormalization group 
equations. Therefore, the relevance of the terms (x) 
and will strongly depend on the bare values 

of their couplin^i^. The latter, being proportional to 
the applied magnetic field, should be controllable in 
experiments^. 


Renormalization group arguments provide an upper 
threshold value on the Luttinger parameter gc = 3/7^ in 
order for a given resonance to be relevant. Estimating 
the actual strength of electron interactions is however 
beyond the scope of the model. More insight may 
be drawn from microscopical approaches such as the 
Hubbard model. This would not only allow to confirm 
the predictions made here concerning the fixed point of 
the theory, but also in principle to observe the transition 
towards the crystal of fractional quasi-particles. In a 
recent papeS^ it has been shown that for /s 
relevant at n > 1, on-site repulsion (however strong) 
is not sufficient in the Hubbard modeP^. This seems 
to imply that a non-zero range interaction is needed 
to observe the predicted effects. Indeed, although the 
bare Hamiltonian of the dot in Eq. includes only 
zero-range interactions, the correction terms pw{x) and 
pz{x) essentially capture the effects of such long-range 
interactions. 

When interactions are not strong enough, the mass of 
r]+{x) does not flow to very large values and additional 
contributions to the density which are negligible in the 
regime treated here, may compete or indeed become 
relevant. In particular, we can expect that lowering 
interactions the Friedel oscillations of fractional quasi¬ 
particles, with (2n -I- 1)Nq/2 peaks, may appear in the 
density. For even lower interactions electrons would 
again become relevant and the standard Wigner and 
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Friedel oscillations, respectively with Nq and No/2 
peaks, would set in. 


IV. CONCLUSIONS 

We have studied the density oscillations of a quantum 
dot created in a spin-orbit coupled quantum wire in the 
presence of a magnetic field and strong repulsive electron 
interactions. When the density is tuned at fractional 
fillings V = kp/kso = l/(2n -I- 1), for sufficiently strong 
interactions a fractional gapped phase develops, for 
which we have shown that the density oscillates at 
the peculiar wavevector 4(2n -I- 1)A:f and thus shows 
(2n -I- l)iVo peaks for Nq electrons in the dot. We 
interpret this strong-interaction result as the formation 
of a crystal of fractional quasi-particles in the ground 
state of the dot. 

We believe that such fractional density oscillations 
can be detected in linear transport experiments on 
state-of-the-art samples, probing the dot with a charged 
AFM tifPUMl. 
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where 


Q'{x) = Q{x) + a-(3+f{x) (A3) 


also includes the effects of commutations between fields. 
Let us first consider the massive sector p = In the 
scaling limit of very large mass, r]+{x) is pinned to the 
minima of the sine-Gordon potential, see Eq. (19). All 


these minima are degenerate, with the j — th minimum 
having energy = — |eo|, Vj with 


A 

So = -— J dx cos[4/(x)] . 


The average of exp [ia+r]^(x)] is then 


(A4) 


/ ^ y 

J—>-oo 


Z.J gt 


^+3 


3^-J 


(A5) 


where po = '^ct+/2, Zj = with a fictitious 

small temperature oc j3~^ which can be set to zero in the 
end of the calculation. If a+ = 21, with I G Z, the above 
expression reduces to (—1)^ If ^ 21 however, the 
partial summations cancel out and the average is zero. 
Thus we conclude that 


OO 

1 — — 00 
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Appendix A: Quantum average of the electron 
density 


In this Appendix we describe the evaluation of the 
quantum average of the electron density p{x) - see 
Eq. (20) on the ground state in the fractional gapped 
phase. 


The long-wave Pun{x) in Eq. (21) term simply evaluates 

No/t 


Each other term - see 
the form 


Eqns. ( |22|24|26| ) 


- can be cast into 


O = cos 


Q{x) + ^ [apPp{x) + / 3 pXp(a;)] 

P^± 


(Al) 


where Q{x) = qx + CQf(x), with q the wavevector of the 
oscillating term and eg a scalar coefficient. To evaluate 
its average, we convert it in the product form 

O = 6**3'(“) ^iPpXpA) + h.c., (A2) 

p=± p—:t 


Note that this sharp condition is due t o the degeneracy 
of the minima of the potential (see Sec. IIIB). 
Furthermore, we have 


(e*^+^+(")) = [E{x)y 


(A7) 


where 

E{x) = ^ 0 (A8) 

due to the large fluctuations of the field x+(x). Thus, 
every operator which contains the field x+(x) will be 
vanishing in the infinite mass limit. 


In the massless sector the (zero temperature) averages 
{r]‘L{x)) and {x-{x)) can be evaluated from the correla¬ 
tors 

Gix,x'-,T) = {Tr'n-ix,T)q-{x',0)) , (A9) 

G{x,x';t) = {TrX-{x,T)x-{x',0)) , (AlO) 

where t is the imaginary time, as 

{ri-{x)) = lim G{x,x;t) (All) 

T-J-a/TTD,, 

iX-ix)) = lim G{x,x-,t) (A12) 

T—^aj'KVj^ 

with Q!(7ru^)“^ a short imaginary time cutoff. Let us 
focus on G{x,x'-,t), as the procedure is analogous for 
G{x,x';t). Introducing the Fourier transform 


Guj{x,x') 


dr 

27r 


G{x,x';T)e 


— lOJT 


(A13) 


0 












from the Hamiltonian in Eq. (|15|) one has 

4Vr, 


9ri VriQri 


Q^{x,x') = 5{x - x'). (A14) 


Imposing ri-{—x) = —? 7 _(x) and rj-(x + 2£) = 77 -(x) as 
appropriate for open-boundary conditions, one finds 


rr ( /^ '^r)9r) ^ sin{qjx) sm{qjx') 

yui(x,x ) — ^ ^ 


C.2+^2^2 


with qj = {Trj£)j. Thus, we get 


{vUx)) = -flog[K{x)], 


where 


K{x) = 


sinh (If) 


^sinh 2 (|f)+sin^(f^)' 
Similarly, one finds 

{X-{x)} = --^log[G(a;)] , 


497 } 


with 


(A15) 


(A16) 


(A17) 


(A18) 


with I an integer. 

Table im shows the coefficients and the wavevector for the 
different oscillating components of the electron density in 

Eq. @. 

It can be seen that p^^\x) = p^^l{x) = — 

(^) = (*) = 0 due to the exponential suppres¬ 

sion induced by the large fluctuations of x+i^)- 
The Wigner harmonics only survive if m = l^n (with 
I S Z*), see Tab. |n] and Eq. (A 6 I, and have wavevector 

b47„A:F. It scales as [if (x)]^ Among them, the most 
relevant term is I = 1 , with wavevector Qo = 47n,fcF- 


Term 


(m 

Pf 


(x) 


7^) 


pT7^) 




77W) 


a+ 


-2 — 

_ In 

77i-7m + l 


-2 — 

_ 7t7, 

1 7m ~7n 


7n—7m + l 


l-7rr 


7n-fm-1 

7n 

7n -7m 

-2ri 


In 

In+lrr 




1-I'n—7m 1 —7m+7Ti 




1 - in-lm 

7n 

7Tr+7m 

Til 




l-7m+7n 

7ti 

Tri^lTm 

7t7 


7n+7m-l In-lm-l 


-2Ln 


-2Ln 


P+ 


-1 


0 


/3- 


-1 


-1 


0 


Wavevector 


2mk^ 


Arnkp 

2(7 

m + Xri- l)fcF 


2(7 

m + Xn- l)fcF 


2(7m-l -f Xn)kF 
2 ( 7 ™ -I- Xn)kF 


2‘iXm 7^ l)i'F 


2 (7m Xn l)i^F 
2 (7m Xn)kF 


2(7m —1 7^)^F 


TABLE II. Coefficients for the different contributions to the 
electron density. 


G{x) = 8 e 7 sinh^ (7j') ^ ^(^) ■ (A19) 


V 2£J 


\2iJ 


Thus in the massless sector one has 
„2 


exp 


exp 


a; 


~{xU7) 


= [if(x)]“-®''/® , (A20) 
= [G{x)f-^^^^^^ . (A2I) 


Summarizing, one obtains for the operator in Eq. (A2) 

c^2 ^2 

(O) = Y.i-iyd^^.21 [Kix)]^ [G(x)]^ 


(A22) 


X [E{x)] ^ cos[(5'(a;)] . 


The scaling limit imposes severe restrictions to Eq. (A22), 
due to the simultaneous requirements a+ = 21, /3_|_ = 0 , 


Going through Tab.[lT|one verifies that similar arguments 
apply for every non-vanishing oscillating term p^j (x): 
they all have the wavevector of the form I ■ dy^/cF (with 

I G Z*) and scale as [K{x)f Thus, we can conclude 
that the only contributions to the electron density oscil¬ 
late with wavevectors I ■ dy^^F- 

Among these terms, the most relevant ones have wavevec¬ 
tor Qo = 47 „A:f and scale as [if(x)]®''^^. When all the 
most relevant terms are collected and physical constraints 
(open boundary conditions and normalization) are im¬ 
posed on {p{x)) to determine Aw and Xz, we obtain 

{ £il 5 'x 

1 - ^ cos [Qoa; - Fi{x)] | , 

(A23) 

where Fi{x) = F 2 {x) = F 3 {x) = 0, ^ 4 ( 0 ;) = 

4{ln - i)f{x), and F 5 (x) = 4 ( 7 „ -b l)f{x). 
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